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<N \ Abstract 
>>'. 

Elastic networks of beams can deform either by stretching or bending of their elements. Rigidity 
of such networks crucially depends on the dominant mechanism of deformation: under the same 
loading conditions, stretch-dominated structures are much stiffer than bending-dominated struc- 
tures. We prove the existence of stiffest networks which deform by extension (or compression) of 
5a \ their elements only, under any loading conditions. Indeed, we provide a convenient set of necessary 

and sufficient geometrical conditions to identify such networks. In addition, values of the effective 

£' 

\^ • elastic moduli of these structures are derived. Our results shed light on the relationship between 

microscopic geometry and macroscopic mechanics of elastic networks, and may be of interest for 
structural applications. 

PACS numbers: 81.05.Zx, 62.20.de, 87.10.Pq, 82.70.Rr 
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Various elastic systems can be understood as networks of interconnected struts which 
deform by a combination of bending, stretching, twisting and shearing mechanisms. Exam- 



ples include polymer gels, protein networ 
atomic lattices and granular materials 7 



cs and cytoskeletal structures 
8|, or foams, wood, and bones 




5], crystal 



Moreover, the pairwise interaction potentials used in standard elastic percolation models 
can also be identified with the strain energy of elastic struts Despite extensive 

research 3), [4, [5, (|, the connection between the mechanical properties of such networks 
on a macroscopic level and the description of their structures on a microscopic level has 
not been completely elucidated yet. Interestingly, under identical loading conditions, some 
structures appear to deform primarily through the local stretching of the struts, while in 
other structures the elastic energy is stored via local bending {^], [nj] (twisting and shearing 
contributions usually are neglected). A transition between these two regimes has been re- 
3orted, both experimentally and numerically, when the mean node connectivity is increased 

In most cases, the apparent stiffness of an elastic network is also dramatically affected 
by the loading conditions: for instance, the square lattice is bending-dominated under shear 
stress, but is stretch-dominated when loaded uniaxially along one of its principal axis. On 
dimensional grounds 1(3], the strain energy density associated with a stretch-dominated 
deformation varies linearly with the relative density 0, defined as the volume of struts per 
unit volume of material, while it scales as <j) 2 for a bending-dominated deformation. Thus, for 
the low-density materials considered here (0 1), a structure deforming primarily through 
the stretching of the beams is much stiffer, and hence much more efficient for structural 
applications. 

In this letter, we analyze the structural conditions under which a network of beams 
deforms exclusively by extension/compression of its elements, for any loading conditions. 
We show that such structures do exist, and we provide a convenient set of criteria on the 
local geometry to identify them. For this purpose, we proceed in two steps: we first establish 
bounds on the elastic moduli of any network of beams, using a variational approach. We then 
show that these bounds are reached for some specific network geometries, determined by a set 
of necessary and sufficient conditions. As expected, these networks deform by axial extension 
or compression of the beams only. We will limit our study to isotropic structures, though 
the reasoning can be transposed without difficulty to materials with lower symmetries. 
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For small strains, the strain energy per unit volume of an isotropic body can be written 
in general as [l5| : 

l kl 

where A is the Lame's first parameter, // the shear modulus (or Lame's second parameter), 
and Um = \ {j^ + f^J are the components of the strain tensor (ut are the components of 
the displacement field u(r)). The elastic properties of any isotropic body are then deter- 
mined by two independent elastic moduli. In what follows we will use the Young's modulus 
E rather than A, which measures the stiffness of a body under uniaxial loading. This elastic 
modulus is defined as E = 4/i (A + fi) / (A + 2fi) and E = fi (3A + 2fi) / (A + fi) for two- 
dimensional (2d) and three-dimensional (3d) materials, respectively. According to Eq. [TJ 
the strain energy is an increasing function of E and \i. 

We first establish bounds on these two elastic moduli using the Principle of Minimum 
Potential Energy (PMPE): consider a body with volume V and prescribed displacements 
on its boundaries; the PMPE states that, among all kinematically admissible displacement 
fields (i.e. all continuous displacement fields satisfying the displacement constraints on the 
boundary), the actual displacement (i.e. the one satisfying the equations of mechanical 
equilibrium) is the one that makes the potential energy S = j v edV an absolute minimum. 

Suppose first that a uniform shear strain is applied on the boundary of an homogeneous 
isotropic body in the xy plan. The macroscopic displacement field associated with such a 
strain is u(r) = (5X/Y)y e x , where X and Y are the dimensions of the sample in the x 
and y directions, respectively, and SX the maximal displacement in the x direction. The 
corresponding strain energy is £ = fi(5X/Y) 2 V/2. On a microscopic level, the body is 
comprised of beams that can have natural curvatures and non-uniform cross-sections. Let 
us denote the beam linking the nodes i and j. We note Sij(l), (/), and iy,(Z) the 
local cross-sectional area, the local tangent unit vector and the position vector originated 
from node i, respectively, where I refers to the curvilinear coordinate along the beam (see 
Fig. Q]). We also note L%j the length of the line segment linking nodes i and j, the unit 
vector along this line segment, and cos 6^- = e a ■ the director cosine of this line segment 
with the a axis (a G {x, y} for 2d materials, and a G {x,y,z} for 3d materials). Let us 
consider the trial displacement field defined as: Uj = (5X/Y)yie x for every node i (jji is the 
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un-deformed 



compression/extension k tJ 
+ rotation 



FIG. 1: Geometry of a typical beam In the initial configuration, the beam can have a 

natural curvature and a non-uniform cross-section. Under the trial displacement field [21 the beam 
is subjected to a translation Uj, a uniform compression/extension kij and a uniform rotation u>ij. 

coordinate of node i along the ?/ axis), and: 

Uij(l) = u { + kijTijty) + u>ij x r y (Z) (2) 

for the point of the central line of beam (i, j) located at position I, where kij = • 
(uj - u.j) /L^ = (8X/Y) cos 6^- cos 6^- and w# = x (uj - u*) /L^ = —{8X/Y) cos 2 ^-e z . 
Thus, each node follows the macroscopic displacement, while each infinitesimal piece of beam 
is subject to a translation Uj, an axial compression/extension and a rotation ujij. Since 
the rotation vector u>jj is constant, there is no bending or twisting of the beam. Moreover, 
we can complete the definition of the displacement field by assuming that the cross-sections 
remain perpendicular to the central line. Then, the only contribution to the strain energy of 
a piece of beam with infinitesimal length dl is the stretching term (Eq/2) Sij(l)kfjdl, where 
E is the Young's modulus of the beam material. For low-density structures, the node con- 
tribution to the strain energy is negligible. Therefore, the energy of the whole structure 
associated with this trial displacement field is E tria i = E /2(8X/Y) 2 'y^v i j cos 2 Q%- cos 2 6^-, 

where the summation is over all the beams that constitute the network, and Vij is the vol- 
ume of beam By comparing £ with StHah an d according to the PMPE, it finally comes 
that the effective shear modulus is bounded as: 

<0<cos 2 ^cos 2 ^), (3) 

Eq 

where the brackets denote the average defined, for any quantity q^, as: (q) = 

V ijQij/ Z2(i,j) V ij- 

We can use the same argumentation to obtain a bound on the Young's modulus: suppose 
that a uniform elongational strain is applied on the boundary of a homogeneous isotropic 



body along the x axis. The macroscopic displacement field associated with such a strain 
is u(r) = (5X/X) (xe x — vye y ) and u(r) = (8X/X) (xe x — vye y — vze. z ) for 2d and 3d 
materials, respectively, where v is the effective Poisson's ratio of the network. In both cases, 
the strain energy associated with such a displacement field is £ = E /2(5X/X) 2 V. We can 
define a trial displacement field similar to Eq. [2J with here the node displacement defined as 
Uj = (5X/X) (x»e x . — vy^y) for 2d materials, and Uj = (5X/X) {xie x — vy^y — uzie z ) for 
3d materials (where Xi, yi, Zi are the coordinates of node i). We can evaluate the parameters 
kij, u)ij, and verify that the only contribution to the associated strain energy comes from 
the axial extension/compression of every infinitesimal piece of beam. Then, the total strain 
energy £ tr iai associated with such a field can be easily evaluated. Using again the PMPE, 
we finally obtain that the effective Young's modulus is bounded as: 



<0^(cos 2 r-z/cos 2 ^) 2 ^) (2d) 



E 

(4) 



-rr < <P ( (cos 2 9 X - v cos 2 0»-u cos 2 9 Z ) 2 ) (3d) . 

The expressions for the bounds [3] and H] involve average quantities which can be evaluated 
using symmetry arguments. For isotropic structures, the expressions of the bounds must be 
identical for any orientation of the applied strain. That is, for 2d structures, the transforma- 
tion 9 X —> 9 X + tp, 6 y —>■ 9 y — ip must leave the expressions unchanged, for any rotation angle 
if. Similarly for 3d structures, the bound expressions must remain unchanged by rotations 
around the three perpendicular axes. Furthermore, bound expressions must be invariant 
under permutation of the axis labels. After a little algebra, these invariance properties lead 
to a set of relations which can be summarized as: 

/ (cos 2 ^) = i <cos 4 r> = ^ 

(5) 

(cos 9 a cos COS6 17 ) = (J3 7^ 7) 
with a, 0, 7 G {x, y} (resp. {x,y,z}), and d = 2 (resp. d = 3) for 2d (resp. 3d) struc- 



tures 16(. The relations [5] imply in particular that (cos 2 9 13 cos 2 6 17 ) = 1/ (d(d + 2)) and 
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cos 9^ cos 6^ = (with (5 ^ 7) . Accordingly, the inequalities [3] and H] simplify to: 





(6) 



JL < A 

" 15 



S, 



< — (8z/ 2 -4i/ + 3) 



(3d) 



It must be noted that the bounds on E are implicit, since they depend on the Poisson's ratio, 
which in turns is related to E and \x by: v = E/2a — 1 (for both two- and three-dimensional 
bodies). 

We now show that inequalities E] become strict equalities for some specific network geom- 
etry determined by the following set of rules (along with the isotropy conditions [5]) : 

a) All the beams are straight. 

b) All the beams have uniform cross- sectional areas: Sy(Z) = 

c) At every junction i of the 2d (resp. 3d) network, and for all e a , Bp, e 7 £ {e^e^} (resp. 
{e x ,e y ,e z }), the following equality is satisfied: 



where the summation is over all the junctions that are connected to the junction i. This 
condition corresponds to a set of 4 (resp. 10) equations per node for 2d (resp. 3d) structures. 

The demonstration is straightforward: according to the PMPE, each inequality in Eq. O 
become a strict equality if and only if the respective trial displacement field is the (unique) 
displacement field that satisfies the equations of mechanical equilibrium. Inspection of force 
and moment balances along each beam and at each junction leads to the three necessary and 
sufficient conditions stated above. Let us make this precise. Consider a specific beam 
we already mentioned that with the trial displacements we chose, each infinitesimal piece 
of beam is only subjected to an axial deformation. Therefore, the force and local moment 
deriving from the trial displacement fields are F»j (/) = —E$Sij (I) fe^ty (I) and My (/) = 0, 
respectively (F^ (I) and M^- (I) are defined as the force and local torque exerted at position I 
by the i-side of the bar to the j-side of the bar). The force balance equation dFij (I) jdl = 
along the beam implies that the direction of F$j does not change with I, i.e. ty (I) = 
for all I. Hence, one also has (I) = const. Therefore, the force balance along every beam 
of the network leads to condition a) and b). The forces acting on any junction i must 




(7) 
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X/E 


K /E 
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2d 


0/8 


0/3 


0/8 


0/4 


1/3 


3d 


0/15 


0/6 


0/15 


0/9 


1/4 



TABLE I: Elastic moduli of stiffest isotropic structures normalized by the Young's modulus of the 
beam material Eq: (//), Young's modulus (E), bulk modulus (k), Lame's first parameter (A), and 
Poisson's ratio (y). 

balance as well. The force exerted by the straight and uniform beam on this junction 
is E Sijkijeij. Thus, mechanical equilibrium at every junction i implies: Y2j s ijkijeij = 0. 
This relation must hold for any orientation of the strain field. By replacing fcy with its 
expression and using the same invariance argument as for the bound expressions, we finally 
obtain, for shear as for elongational strain, that condition c) is also required. Obviously the 
moment balance equation rjj (I) x F^- (I) = along every beam is also satisfied for 
structures satisfying conditions a) and b). 

Let us discuss some of the implications of the conditions a)-c). First, as long as the 
isotropy conditions (jSJ) are satisfied, the two elastic moduli attain their respective bounds 
under the very same conditions. We can thus eliminate v and express the value of the 
Young's modulus: E = E Q (j)/3 (resp. E = E <f)/G) for 2d (resp. 3d) networks. For isotropic 
structures, any other elastic modulus is related to E and /x, and hence can be easily evaluated. 
Values of most common moduli are reported on table [B 

Interestingly, the Poisson's ratio of such structures is less than the one expected for an 
incompressible body (1 for 2d body, 1/2 for 3d body). It can also be noticed that the 
two Lame's parameters A and /i are equal. Values of the shear, bulk, and Young's moduli 



or low-density (0 <C 1) 



reported on table [I] coincide with the Hashin-Shtrikman bounds [itJ 
2d structures, but are strictly lower (i.e. sharper) for 3d structures [l8j. 

It is worth mentioning that condition c) implies in particular that the vectorial equality 
^2 s ij e ij = is satisfied in stiffest structures. Therefore, these structures have maximal bulk 



modulus 



21 



22 



231 ] as well, for a 



201 ] and maximal electrical (or thermal) conductivity 
given value of 0. 

Our results also reveal why open-cell foams are almost always bending-dominated. Foams 
are particular cellular materials: their structures result from a surface minimization process, 
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leading to geometrical and topological rules, known as Plateau's laws. Assuming that edges 
have identical cross-sections, Plateau's laws state that edges meet in threefold (resp. fourfold) 
junctions with equal angles of 120° (resp. 109.5°) in 2d (resp. 3d) foams. One can easily 
check that such node configurations cannot satisfy condition c) [24j. Therefore, open-cell 
foams will rather deform by bending of their edges. 

Some examples of stiffest structures can be given. Let us inspect the particular case of 
periodic structures. 2d isotropic structures with periodic arrangement must have hexagonal 
symmetry 15j. If we restrict our attention to networks with identical beams, only two 



candidates are left: the hexagonal and triangular lattices. Although the hexagonal lattice 
satisfy the isotropy conditions 0, it does not satisfy the condition c), and therefore is not 



stretch-dominated, as it as been already reported 



191 ] [actually, it can be shown that a 



three-fold junction of beams with equal sections cannot satisfy condition c)]. The triangular 
lattice, on the other hand, satisfies all the conditions. Indeed, the elastic moduli of this 
lattice match the values expected for stiffest structures [l^l, 

The 3d case is less obvious. The structure depicted on Fig. [2^ (where the bars join the 
center of the faces of a regular dodecahedron to its center) is an example of 3d isotropic 
stiffest structure: it satisfies the isotropy conditions [5] and the set of criteria a)-c). However, 
it cannot be repeated periodically because the dodecahedron does not perfectly tile the space. 
In fact, we could not identify any periodic and isotropic rigid structures with one single kind 
of bar. However, we found stiff structures with two types of bars. For instance, the structure 
depicted on [2b satisfies all the required conditions and can be repeated periodically: the 
bars join the centers of the faces of a Kelvin cell (which is known for tiling the space). The 
length ratio of the two kinds of bars is imposed by the Kelvin cell geometry: I2/I1 = VS/2. 
We then adjust the section ratio for the isotropy conditions [5] to be satisfied. We obtain: 
S2/S1 = 3^3/4. 

In summary, we established a convenient set of conditions on the local geometry of 
isotropic networks of beams to be stiff under any loading conditions. The elastic mod- 
uli of these networks are also derived, and can be simply expressed in terms of the Young's 
modulus of the beam material and the relative density of the structure. Examples of two- 
and three-dimensional stiff structures are also given. These results explain why most cellular 
structures, and foams in particular, are bending-dominated. We hope these results are of 
interest for structural applications as well as for our understanding of biological systems. 
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FIG. 2: (color online) stiff three-dimensional structures: (a): structure with identical bars: the 
bars join the center of the faces of a regular dodecahedron to its center. This structure satisfies 
the isotropy conditions [5] and the set of criteria a)-c) but cannot be repeated periodically, (b): 
structure with two types of bars. The bars join the centers of the faces of a Kelvin cell. Such 
structure satisfies both the isotropy conditions and the set of criteria a)-c) and can be repeated 
periodically. 

We thank J. M. Cornelissen, A. Emeriau and C. Gay for critical reading of the manuscript. 
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